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We study the Coulomb drag between two single graphene sheets in intrinsic and extrinsic graphene 
systems with no interlayer tunneling. The general expression for the nonlinear susceptibility appro- 
priate for single-layer graphene systems is derived using the diagrammatic perturbation theory, and 
the corresponding exact zero-temperature expression is obtained analytically. We find that, despite 
the existence of a non-zero conductivity in an intrinsic graphene layer, the Coulomb drag between 
intrinsic graphene layers vanishes at all temperatures. In extrinsic systems, we obtain numerical 
results and an approximate analytical result for the drag resistivity pr>, and find that pd goes as T 2 
at low temperature T, as 1/d 4 for large bilayer separation d and 1/n 3 for high carrier density n. We 
also discuss qualitatively the effect of plasmon-induced enhancement on the Coulomb drag, which 
should occur at a temperature of the order of or higher than the Fermi temperature. 

PACS numbers: 



I. INTRODUCTION 

With the recent advent of the experimental fabrica- 
tion of a single layer of graphene, the electronic and 
transport properties of this newly discovered material 
have been intensively studied both experimentally 1 -^ and 
theoretically^. Whereas electronic structure experi- 
ments have revealed detailed subtle many-body effects 
on the graphene energy spectrum, transport experiments 
have also revealed some apparently unusual features of 
graphene transport properties, most noticably, that the 
conductivity has a non-zero minimum value around zero 
bias gate voltage. Up to now, the transport experiments 
performed have been focused only on the longitudinal 
and Hall transport properties (in both weak and strong 
magnetic fields, including the quantum Hall regime) and 
weak localization, where all of these phenomena depend 
essentially only on the physics of scattering of individ- 
ual single quasiparticle from impurities with electron- 
electron many-body interaction effect playing the role of 
a small quantitative correction. In two-dimensional elec- 
tron gas (2DEG) semiconductor bilayer structures (e.g. 
modulation-doped GaAs/Al^Gai-^As double quantum 
wells), electron-electron scattering between the 2DEG 
layers give rise to the effect of Coulomb drag where a 
'drag' current is induced purely from the momentum ex- 
changes through interlayer electron-electron scattering 
events. One measures the effect of Coulomb drag by 
the drag resistivity pu, which is defined by the induced 
drag electric field in the open-circuited passive layer per 
unit applied current density in the active layer. In high- 
mobility samples where the disorder is weak, pd goes as 
T 2 at low temperatures T, and as 1/d 4 for large bilayer 
separation d (Refs. yla). 

In this paper, we investigate the Coulomb drag in 
graphene "bilayer" systems with no interlayer tunneling, 
considering both the intrinsic (chemical potential fi = 0) 
and extrinsic (fj, ^ 0) cases. We emphasize right in the 
beginning so that there is no semantic confusion what 



we mean by the terminology 'bilayer' graphene. Our 
'bilayer' graphene is two isolated parallel 2D graphene 
monolayers separated by a distance d, with no interlayer 
tunneling. The electronic structure of each graphene 
monolayer is thus unaffected by having the other layer. 
Each graphene layer is assumed to have its own vari- 
able carrier density in the extrinsic case. Our system is 
thus different from the ordinary bilayer graphene where 
d ~ 1 — 5 A with strong interlayer tunneling. Throughout 
this paper, we shall also use the terms "undoped" and 
"doped" interchangeably with "intrinsic" and "extrin- 
sic" respectively; keeping in mind that in experiments 
the chemical potential can be changed by both chem- 
ical doping and gating with an applied voltage. The 
Coulomb drag in graphene is interesting not only be- 
cause it is a novel material with a linear energy spec- 
trum, but also because it only spans a thickness of a 
single carbon atom, the electrons are much more con- 
fined along the perpendicular direction comparing with 
2DEG in a quantum well, where the finite width thick- 
ness has to be taken into account in any quantitative 
comparison with experiments. Thus, the Coulomb drag 
phenomenon in graphene is expected to be theoretically 
very well accounted for with two zero-thickness graphene 
sheets. In addition, tunneling is only appreciable when 
the the out-of-plane it orbitals from the two graphene 
sheets start to overlap with each other at an interlayer 
distance d of a few angstroms (d ~ 3. 5 A in naturally 
occuring graphite), making it possible to study the ef- 
fect of Coulomb drag for an interlayer separation d down 
to a few tens of angstroms, about an order of magnitude 
smaller than is possible in the usual double quantum well 
systems without appreciable tunneling. 



II. FORMALISM 

Graphene has a real-space honeycomb lattice structure 
with two interpenetrating sublattices, giving a honey- 
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FIG. 1: Diagrams contributing to the drag resistivity Eq. 0. 
The double wavy lines represent the screened interlayer 
Coulomb potential Eq. §4$ and the vertices on the left and 
on the right denote charge current in the two layers. Dashed 
vertical lines next to the vertices denote impurity vertex cor- 
rection to the charge current. 



comb reciprocal space structure with the Brillouin zone 
forming a hexagon. At each of these corners where the 
high symmetry point is denoted as K or K', the energy 
dispersion is linear and the effective low-energy Hamilto- 
nian at the vicinity of these K and K' points is given by 
Hq = v& ■ fc. The eigenenergy is e k \ = Xvk, where the 
chirality label A = ±1 describes conduction band (A = 1) 
and valence band (A = —1) spectra. Throughout this 
paper we choose h = 1 unless it is written out explicitly. 
The transformation which diagonalizes the Hamiltonian 
is a local transformation, dependent on the momentum 
fc through tan^fe = k y /k x : 
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The central quantity in the Coulomb drag problem is 
the nonlinear susceptibility 7 '^, T: 
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{ [Gt q (e) - G%_ q (e)] G£(e + ui) J(k)G%(s + u>)} 
+ {q,u^-q,-u}, (2) 

where G^' (e) = (e — Hq ± i/2t)^ 1 denotes, within the 
Born approximation for the self-energy, the impurity- 
averaged retarded/advanced Green function, t the life- 
time due to impurity scattering, J the charge current ver- 
tex, np the Fermi function and 'tr' the trace. In the rest 
of this paper, we shall simply denote the x-components 
of r and J as F and J. The drag conductivity, diagram- 
matically shown in Fig. [TJ 
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here subscripts T' and '2' are the labels for the two single- 
layer graphene, U\ 2 is the screened interlayer potential, 
which in the random phase approximation (RPA) is given 
by 



U12 (q,u) 



(4) 
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where d is the interlayer spacing, V(q) = 2-Ke 2 /q is 
the bare Coulomb potential, n (q, ui) is the graphene 
polarizability^. Instead of the drag conductivity, in ex- 
periments one measures the drag resistivity pr> which is 
defined as (here W and L are the width and length of 
the sample, respectively) pr> = (W/L)(V 2 /Ii), the ra- 
tio of the induced voltage in the passive layer V 2 to the 
applied current in the active layer I\. The drag resis- 
tivity is then obtained from the drag conductivity as 
Pd = 00/(0X10x2 - 0d) ^ 00/0X10x2, where cr L i,2 is 
the longitudinal conductivity of the individual layer 1 or 
2. 

In this paper, we shall restrict ourselves to the Boltz- 
mann regime {lot > 1 or q( > 1, where I = vt is the 
mean free path) corresponding to weak impurity scat- 
tering, which is the case relevant to actual experimen- 
tal situations where high-mobility samples with dilute 
impurities are used. The longitudinal current for the 
graphene Hamiltonian is J — edHo/dk = ever. In the 
presence of impurities, vertex correction to the current 
is taken into account within the impurity ladder approx- 
imation, which gives the impurity-dressed current ver- 
tex as J — (rt r /r)ewo - , where the transport time Tt r for 
graphene is given by 



: : = ■K^2n i \u l {k-k')\ 2 {l-cos 2 6 + )5{e k+ -e k i + ) 



k—kj. 



with m and Ui being respectively the impurity density 
and impurity potential; 8+ = 4>k+q — fik being the scat- 
tering angle from momentum fc to fc + q. We have given 
Eq. ([5]) for the transport time here only for the purpose of 
completeness, as our final results for the drag resistivity 
do not depend on r tr . 

Following Ref. i, we express the nonlinear susceptibil- 
ity Eq. © as 

r{q,u) = T ]T ^[j AA (fc + g)- Jw(fc)] (6) 



A,A'=± fc 



Im\ (1 + AA'cos6»+) 



CfcA' — Cfc+qA 



where J(k) — U^JUk is the impurity-dressed charge 
current vertex expressed in the chiral basis. Eq. ([6|) 
is different from the nonlinear susceptibility for regular 
2DEG with quadratic spectrum in two ways: (I) there 
are contributions to the electron-hole excitations com- 
ing from intraband transitions (A = A') and interband 
transitions (A 7^ A'); (2) there is an additional factor 
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(1 ± cos£> + )/2, which derives from the Berry phase struc- 
ture of the graphene Hamiltonian. Furthermore, the 
nonlinear susceptibility Eq. (|6|) is not directly propor- 
tional to the imaginary part of the polarizability as in 
regular 2DEG, because here the current J(k) is not di- 
rectly proportional to the momentum k. Eq. ((6J has 
the same formal structure as in the case of a regu- 
lar 2DEG with Rashba/Dresselhaus spin-orbit coupling 9 , 
where A = ±1 describes the two spin-split bands. The 
finite off-diagonal components of J do not contribute to 
the nonlinear susceptibility, and only the diagonal com- 
ponents J\\ — A(r tr /r)ez; cos^fe enter into the expression 
Eq. ([6]), corresponding to electrons moving in the conduc- 
tion band (A = 1) with a velocity of constant magnitude 
{j\ z /t)v and valence band (A = —1) with — {t^/t)v. In 
the following, we consider the Coulomb drag between in- 
trinsic graphene layers and extrinsic graphene layers sep- 
arately. 



III. DRAG IN INTRINSIC GRAPHENE 
SYSTEMS 



lows. We first make a change of the integration variable 
k' = —k in Eq. ([S]), and use time- reversal symmetry to 
obtain 

T(q,u)=T £ E [-Jxx(k'-q)+Jyy(k') 



A,A'=± k' 



Im\ (l + AA'cos0_) 



n F{e k >y)-n F {e k ,_ qX ) 
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where 9- = <pk — (j>k—q- This is so far general. Next we 
impose the symmetry requirements of the bands about 
e = 0, i.e. e fe / )A = -e h >_ x and J k , x = -J fe ',_ A , and 
then change the band labels as r' — —A, r = —A' in 
Eq. 0. Finally, using the relation np(e k ' _ r ) = 1 — 
n ^( e fc'.r) va hd for the undoped case /i = 0, we arrive 
at r = — r, i.e. T(q, u) = 0. This result holds true for 
any type of spectrum where the two bands have a mirror 
symmetry across e = 0, and any bilayer system with one 
or both of the layers having such a band symmetry with 
zero doping always results in an overall vanishing drag at 
all temperatures. 



For the case where the graphene layers are undoped 
fx = 0, we first state the main result: the drag conduc- 
tivity between two intrinsic graphene layers, or between 
one extrinsic and one intrinsic graphene layers, is identi- 
cally zero. This is not at first sight a trivial consequence 
of zero doping if one recalls there is a finite conductivity 
(so-called the "minimum conductivity") at zero doping 
in graphene. A physical explanation and a general ar- 
gument for the reason why this is so is in order. When 
the Fermi level is at the Dirac point, the only process for 
electron-hole pair creation will be interband electron ex- 
citation from the valence band to the conduction band by 
which equal numbers of electrons and holes are created. 
In the mechanism of Coulomb drag, the applied electric 
field drives the electrons (or holes) in the active layer in 
the, say, positive (negative) direction; through Coulomb 
scattering, momentum is transferred to the passive layer, 
which drives the carriers (regardless of whether these are 
electrons or holes) in the same direction as the momen- 
tum transfer. In doped systems where there is only one 
type of carrier (either electron or hole), this gives a fi- 
nite drag current in the passive layer. Now, in undoped 
systems where a perfect electron-hole symmetry exists, 
there are two cases for consideration: (1) If the active 
layer is undoped, equal numbers of electrons and holes in 
the active layer will be driven in the opposite direction by 
the applied electric field, and the net momentum transfer 
is thus zero. There will be no drag regardless of what the 
passive layer is. (2) If the active layer is doped while the 
passive layer is undoped, equal numbers of electrons and 
holes in the passive layer will be driven in the same di- 
rection by the momentum transfer, therefore resulting in 
a vanishing drag current. The conclusion of these qual- 
itative considerations amounts to a vanishing nonlinear 
susceptibility T(q,Lu) — 0, which we now proof as fol- 



IV. DRAG IN EXTRINSIC GRAPHENE 
SYSTEMS 

We now move on to the drag between finite-doped 
graphene layers. We first provide the exact analytical re- 
sults for the nonlinear susceptibility Eq. ^ with Fermi 
energy ep > (in the following x — q/kp and y — uj/ep, 
and f = r/(2efc F r/vr)): 



f (x, y) = f in tr a 0(a: - y) + f intcl 9(y - x), 



(8) 



where 9 is the unit step function, ri ntra is the intraband 
contribution to the nonlinear susceptibility given by the 



terms with A = A' 



Eq. ©, 



ftotrafotf) = ^\/z 2 - V 2 {2y/(y + x- 2)(y - x - 2) 



- \J x 1 - y 2 



tan 



y/(y + x - 2)(y - x - 2) x /x 2 - 
x 2 - 2 - (y - 2)y 



+7T0 [y(y -2)~x 2 + 2]]}9(2-x-y)-{y^-y}, (9) 

and r; ntcr is the interband contribution given by the 
terms with A ^= A' in Eq. ((6|), 



^y 2 - x 2 [2^(x + y-2)(x-y + 2) 



tan 



^(x + y-2)(x-y + 2) x /y 2 -x 2 
x 2 - 2- (y- 2)y 



-ltd [x 2 - 2 - (y - 2)y]]} 9{x + y - 2)9{x -y + 2). (10) 

The intraband contribution correponds to electron-hole 
excitations in the vicinity of the Fermi level within the 
conduction band, which occur at uj < vq; whereas the 
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interband contribution corresponds to electron-hole ex- 
citations from the valence band to the conduction band, 
which occur at uj > vq. Using Eqs. ((3]), (|9|)- (fT0| and 
the expression for the graphene polarizability^, we have 
calculated numerically the drag resistivity pp> for differ- 
ent values of interlayer distance d and density n (Fig. [3]). 
Before we proceed to discuss our numerical results, it is 
instructive to obtain an analytical formula for the drag 
resistivity under certain approximations. To this end, we 
first define the standard dimensionless parameter for the 
excitation energy u — y/x = uj/vq in the Fermi liquid 
theory, which is the ratio of the phase velocity of the 
excitation to/q to the quasiparticle velocity v. At low 
temperatures and with large interlayer separation, the 
dominant contribution to the drag conductivity Eq. Q 
comes from region with small q and to, consequently the 
nonlinear susceptibility can be evaluated in the limits of 
long wavelength x <C 1 and low energy m«1, allowing a 
closed-form expression for T(q, ui) to be extracted. The 
interband (A ^ A') contribution in Eq. ^ is in general 
smaller than the intraband (A = A') contribution by an 
order 0(x 2 ), and vanishes in the limit u -C 1 as seen from 
Eq. ((SJ. This is because, in the presence of a finite Fermi 
level, electrons take more energy to transition from the 
valence band to the conduction band (interband) than to 
transition within the conduction band (intraband), and 
with a small excitation energy the channel of interband 
transition becomes inaccessible. With the above assump- 
tions, Eq. can be evaluated as 



r(x, y = ux) = r intra = -Au 



- (1 + ux) 



t-9(l-\t-\) 



1-t 2 



(1 — ux) 



-4ux 



M(i-IM) 
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where t± = u± x(l — u 2 )/2. Eq. (fTl"]) is larger than the 
corresponding expression for the nonlinear susceptibility 
in regular 2DEG by a factor of 4, due to an extra 2x2 
degrees of freedom coming from the spin and valley de- 
generacies in graphene, in addition to the two sublattice 
degrees of freedom which give rise to the conduction and 
valence bands. 

The longitudinal conductivity can be obtained from 
the impurity-dressed current J = (T tr /r) ev& using the 
Kubo formula to give ox = e 2 vD 1 where v — 2kp/irv is 
the graphene density of states and D = v 2 T tI /2 is the 
diffusion constant. This Kubo formula result is identical 
with the Boltzmann theory result ox = (e /h)2spTt T /h. 
Incidentally, for short-range impurities the transport 
time is related to the lifetime simply by r tr = 2r due 
to the suppression of backscattering from impurities in 
graphene. There are two types of disorder in substrate- 
mounted graphene, one being the charged impurities 
coming from the substrate; and the other being the neu- 
tral impurities intrinsic to the graphene layer itself. In 
theory, the type of disorder essentially boils down to the 



expression of the transport time r tr in ox, which yield 
different types of functional dependence of the conduc- 
tivity <xl on the carrier density. In experiments, the con- 
ductivity is observed to increase linearly with density, a 
fact alluding to the dominance of the charged impuri- 
tity scattering in substrate- mounted graphene samples. 
We emphasize that this dependence on different types of 
disorder does not affect the expression of the drag re- 
sistivity as the transport time Tt r is explicitly canceled 
out between ox> and 0x10x2- Therefore, our calculation 
and conclusions apply equally to bilayer systems with 
substrate-mounted (where charged impurity scattering 
plays the more dominant role) or suspended graphene 
samples (where there is only neutral impurity scatter- 
ing). 

In the expression of the drag conductivity Eq. (J3J), the 
dominant contribution of the integral comes from the re- 
gion where qd < 1, and for large interlayer separation d 
satisfying d^ 1 <C k F , 9tf , the interlayer potential Eq. ^ 
can be approximated as U12 — g/47re 2 sinh(<7<i)riin2. 
Furthermore, the denominator sinh 2 (lo /2k bT) in Eq. ([5]) 
also restricts the upper limit of the u> integral to a few 
~ fcsT, therefore at low temperatures only small values 
of w/ej? contribute to Eq. (J3j) . As a consequence, the 
polarizability for doped graphene can be approximated 
by the static screening result H(q, u>) — v. Now, using 
Eq. (fn) for the nonlinear susceptibility, the drag resis- 
tivity is obtained as 



PD 



h ttC (3) (k B T) 2 
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e 2 32 EFieF2 {qTFid)(qTF2d) (k F id)(k F 2d)' 

where qtf = 4e 2 fc^ jv is the Thomas- Fermi wavenumber 
for extrinsic graphene^. The drag resistivity Eq. (112p . 
valid for low temperatures T <C T F and high density 
and/or large interlayer separation k F d 1, has exactly 
the same form as in the regular 2DEG drag, exhibiting 
the same dependences of temperature (~ T 2 ), interlayer 
separation (~ l/d 4 ) and density (~ (niri2)~ 3 / 2 )- 

Our numerical calculations and analytical results 
Eq. (I12p are compared in Fig. [51 showing that Eq. (TT2")) 
becomes an increasingly accurate approximation to the 
full numerical results with increasing values of kpd. The 
fact that the exact numerical results shown in Fig. [2] dis- 
agree more strongly with the analytic result of Eq. (TT2"]) 
for smaller values of k F d is understandable, since the 
analytic formula given in Eq. (|I2p applies only in the 
asymptotic kpd 3> 1 limit, and for lower carrier den- 
sity and/or interlayer separation, Eq. (|I2p simply does 
not apply. In particular, for kpd = 1, the exact nu- 
merical result for Coulomb drag is a factor of 4 larger 
than that given by Eq. (fT2]) . This trend of an increas- 
ing quantitative failure of the asymptotic analytic drag 
formula for lower values of kpd has also been noted in 
the literature^ in the context of low-density hole drag in 
bilayer p-GaAs 2D systems. For small kpd, backscatter- 
ing effects in Coulomb drag, which are unimportant for 
kpd 3> 1, become important. 

On the other hand, our numerical results also show 
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FIG. 2: (Color online) p D /(T/T F ) 2 as a function of T/T F 
for Coulomb drag between two identical extrinsic graphene 
sheets, with values of k F d = 10 (solid lines), 5 (dashed lines), 
and 1 (dot-dashed lines). Numerical results are indicated 
with bold (black) lines and analytical results Eq. (I12|l with 
thin (grey/red) lines. The analytical results become an in- 
creasingly accurate approximation to the full numerical re- 
sults with increasing k F d (i.e. increasing n or d). 



that the temperature dependence of pu remains very 
close to T 2 within a wide range of temperatures for 
typical experimental values of d and n (e.g. kpd = 5 
with n = 5 x 10 n cm~ 2 and d ~ 400A). The ratio of 
the Fermi temperature for graphene to that for regular 
2DEG with parabolic spectrum (with effective mass m) is 
Tp (graphene) /Tp(2DEG) = m»/fiy , ro J so for low den- 
sities e.g. n = 10 n cm -2 , Tp = 430K for graphene can 
be larger by an order of magnitude than Tp = 42K for 
GaAs 2DEG, and the temperature dependence of pd for 
graphene drag therefore remains very closely T 2 up to 
about several tens of kelvin where the low temperature 
regime T <C Tp still remains valid, whereas for drag in 
regular 2DEG systems departure from the T 2 dependence 
of pd typically occurs at T < 10K. The drag resistiv- 
ity is calculated numerically for various values of d and 
n and higher values of temperature up to T = 0.2Tp 
(Fig. [3]); /?D is seen to grow slower and slower than T 2 
as temperature increases. Similar dependence on tem- 
perature is also observed for drag in regular 2DEG bi- 
layer systems before T reaches > 0.2Tp, beyond which 
Pd/T 2 starts to increase due to plasmon enhancement to 
the drag resistivity^. We discuss the effect of plasmon 
enhancement to the Coulomb drag in graphene bilayer 
systems in the following. 

In regular 2D bilayer systems, enhancement to the drag 
resistivity due to coupled plasmon modes comes into play 
with increasing temperature^. There exist two plasmon 
modes, the so-called optic and acoustic modes, for which 
the electrons on the two layers move collectively in phase 
and out of phase, respectively, with each other. The en- 
ergy dispersion lines for these plasmon modes lie above 
the electron-hole excitation continuum (i.e., the region of 
uj vs. q where the imaginary part of the polarizability is 
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FIG. 3: p /(T/T F ) 2 vs. T/T F for higher values of T up to 
0.2Tf. Upper panel: for fixed interlayer distance d — 500 A 
and different values of density n = 10 11 cm -2 (solid line), 
5 x 10 11 cm -2 (dashed line), 10 12 cm -2 (dot-dashed line), cor- 
responding to T F ~ 431K, 963K, 1361K respectively; lower 
panel: for fixed density n = 10 11 cm -2 and different values of 
interlayer distance d = 300A (solid line), 150A (dashed line) 
and 30A (dot-dashed line). 



non-zero, ImTl(q, u>) ^ 0) at zero temperature, and are 
not excited at low temperatures. They can be excited, 
however, at higher temperatures when the electron-hole 
excitation continuum occupies higher values of the exci- 
tation energy to, a consequence of the increasing gradi- 
ent with increasing momentum k in the parabolic energy 
dispersion relation. Absorption or emission of a plas- 
mon can occur when the electron-hole excitation con- 
tinuum starts to overlap with the plasmon dispersion. 
On the other hand for graphene, because the gradient 
of the linear dispersion relation is constant, increasing 
temperature does not increase the range of the possible 
intraband excitation energies, the electron- hole excita- 
tion continuum being always bounded by lu < vq. This 
means that the plasmon excitation energy will always 
be out of reach from the intraband excitation channel 
at all temperatures. However, the case is different with 
the interband excitation, for which the electron-hole ex- 
citation continuum overlaps already at T = with the 
plasmon dispersion at about*^ u> > £p. This means that 
plasmon-induced enhancement of the drag resistivity in 
graphene occurs, solely due to interband transitions, at 
a temperature T > Tp; whereas for regular 2D systems 
plasmon-induced enhancement occurs already before T 
reaches Tp (at about T ~ 0.5Tp). 



V. CONCLUSION 

In conclusion, we have formulated the Coulomb drag 
problem for graphene bilayers. The drag resistivity is 
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zero for intrinsic graphene. For extrinsic graphene, the 
interband contribution to the drag due to electron-hole 
excitations is suppressed at low temperatures, and the 
Coulomb drag is due predominantly to the intraband con- 
tribution near the Fermi surface in the conduction band. 
We have obtained exact analytical results at T = for 
both intraband and interband contributions to the non- 
linear susceptibility, and obtained the drag resistivity nu- 
merically. We have also derived an approximate analyt- 
ical result for the drag resistivity valid for low temper- 
atures, high density and/or large interlayer separation. 
We find both similarities and differences for the graphene 
drag resistivity compared with that for regular 2DEG 
with quadratic energy spectrum. At low temperatures, 
graphene drag resistivity exhibits the same temperature, 
bilayer distance and density dependences as regular 2D 
systems. For low densities n < 10 cm , the low tem- 
perature regime where the T 2 dependence of the drag re- 
sisitivity holds extends by an order of magnitude that for 
regular 2D systems, as the Fermi temperature is higher 
for the same carrier density in graphene than in regular 
2D systems. In contrast to regular 2D bilayer systems, 
there is no contribution to plasmon-induced enhance- 
ment of the drag resistivity due to intraband excitations, 
and the only contribution to plasmon-induced enhance- 



ment comes from interband excitations, which occur at 
temperatures T > sp. The coupled plasmon modes in 
graphene bilayer systems can therefore be probed ex- 
perimentally with drag resistivity measurements at high 
enough temperatures or at low densities. 

Finally, we comment on the possible effect of disorder 
which has been argued^ to be important for substrate- 
mounted graphene monolayers in the low carrier density 
regime. In particular, the low-density graphene monolay- 
ers will be dominated^ by spatial inhomogeneities asso- 
ciated with electron-hole puddles induced by the charged 
impurities in the substrate. In the presence of such den- 
sity inhomogeneity, the Coulomb drag in the low-density 
regime may deviate substantially from our theory based 
on the spatially uniform carrier density (in each layer) 
model. In particular, there could be large drag fluctu- 
ations including even negative drag in this low-density 
regime. 
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